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Приведен обзор результатов исследований и новые экспериментальные данные авторов по изменениям изотопного состава в 
приповерхностных слоях твердых тел в процессах ионного распыления, ионного рассеяния, вторичной ионной эмиссии, ионной 
имплантации, электролизного насыщения изотопами водорода, термодиффузии из внешнего источника. Во всех перечислен- 
ных процессах наблюдаются заметные изменения естественных изотопных распространенностей, как в составе приповерхност- 
ных слоев, так и в составе вторичных (эмитированных, отраженных) частиц. Исследования выполнены методом вторичной ион- 
ной масс-спектрометрии и энерго-масс-спектрометрии вторичных ионов на образцах молибдена, никеля имплантированного 
медью, на различных модификациях титана, на тонкопленочных системах титан-алюминий и других. Выявлены общие законо- 
мерности в изменениях изотопного состава в различных процессах. Обсуждаются возможные механизмы этих изменений. 


Введение 

В статье выполнен обзор опубликованных в 
последние годы работ и приведены новые экспери- 
ментальные данные по изменениям изотопного 
состава в приповерхностных слоях твердых тел в 
ходе ионного распыления, отражения ионов от по- 
верхности, вторичной ионной эмиссии, ионной 
имплантации, электролизного насыщения изото- 
пами водорода, термодиффузии из внешнего ис- 
точника. Такие исследования обусловлены потреб- 
ностью материаловедения, энергетики и других 
прикладных отраслей знаний, использующих тех- 
нологии модификации приповерхностных слоев, и 
занимающихся изучением механизмов миграции 
атомов, трансформации изотопного состава при 
различных воздействиях на поверхность. Эти ис- 
следования представляют и академический инте- 
рес как примеры модельных объектов неравновес- 
ных физических систем. 

Изменения изотопного состава в ионных пуч- 
ках и приповерхностных слоях твердых тел изуча- 
лись во многих работах по ионному распылению и 
вторичной ионной эмиссии [1—10] , по рассеянию 
ионов от поверхности [11, 12], при водородном на- 
сыщении и термодиффузии из внешнего источни- 
ка [13-20], при ионной имплантации [21-23]. В 
этих явлениях важно выявить основные процессы, 
обусловливающие изменения изотопного состава в 
приповерхностных слоях твердых тел. 

Экспериментальные результаты и обсуждение 

1. Изотопные эффекты при ионном распылении и во вто- 
ричной ионной эмиссии 

Первыми работами, в которых исследовался 
изотопный эффект (ИЭ) во вторичной ионной 
эмиссии (ВИЭ), являются, по-видимому, работы 
авторов [1], хотя обнаружение предпочтительного 
распыления (ргеіегепііаі $рийегіп§), то есть изотоп- 
ного эффекта в распылении относят к 1977 году 
(см. [2] и библиографию этой работы). В [1] обна- 
ружено, что в ВИЭ ряда минералов, металлов, ок- 


сида меди и сплава ОаАІ ионами кислорода отно- 
шение выхода ионов легкого изотопа к тяжелому в 
энергетическом диапазоне 0...20 эВ выше соответ- 
ствующего стандартного отношение распростра- 
ненностей изотопов. Наблюдается обогащение 
"облака” распыленных ионизованных частиц лег- 
кими изотопами. Степень обогащения зависит от 
атомного номера элемента, от матрицы, в которой 
элемент находится, и скорости ионов. 

Практически одновременно с работами [1,2] по- 
явилась работа [3], где исследовано пространствен- 
ное и энергетическое распределения величины из- 
меряемого в ВИЭ изотопного отношения. Работа [3] 
подтвердила выводы работ [1] относительно обога- 
щенности ВИЭ легкими изотопами и зависимости 
изотопного отношения вторичных ионов одного и 
того же элемента от их энергии. Авторами [3] под- 
робно исследовались для ряда чистых металлов за- 
висимости/ ІЯ (-Е)=Л^ І (^/А + й (^, где № Ь (Е) и Л%( Е) 
- энергетические спектры вторичных ионов легкого 
и тяжелого изотопов. Было установлено, что в зави- 
симости /, я ( Е) в исследованном диапазоне энергий 
(от 0 до 80 эВ) имеется пик или плато при разных 
энергиях для разных металлов. Как правило, этот 
пик расположен в области низких энергий для лег- 
ких элементов и смещен к высоким энергиям для 
тяжелых. Подробно вопрос о зависимости /ін(Е) в 
диапазоне Е= 4.. .30 эВ изучался также в работе [4] с 
целью проверки предложенной в этой работе моде- 
ли эффекта. В дальнейшем закономерности ИЭ в 
распылении и ВИЭ исследовались в [2, 4-7 и др.]. 
Для теоретического объяснения ИЭВИЭ в [1-7] 
сделано следующее. Привлечены известные форму- 
лы для вероятности ионизации (в работе [3]) и экс- 
периментальная зависимость вероятности иониза- 
ции от скорости атома (ѵ) в виде ехр(-ѵ 0 /ѵ), ѵ 0 - 
константа (в работах [2,4]). Эффект изменения изо- 
топного состава получается как следствие зависи- 
мости вероятности ионизации от массы атома. Это, 
однако, не объясняет физической сущности меха- 
низма эффекта. Целью наших исследований [9, 10] 
являлось экспериментальное изучение ИЭ ВИЭ и 
попытка углубления представлений о его механиз- 
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ме. Экспериментально изучены зависимости /, я ( Е) 
для металлов, круг которых частично перекрывается 
с теми, что исследовались в [1, 3-7]. Получено каче- 
ственное подтверждение экспериментальных выво- 
дов указанных работ, хотя в количественном отно- 
шении для ряда металлов у нас наблюдался эффект 
превышающий величины, опубликованные ранее, 
как видно, например, из рис. 1. 

І 92 МоЛ 10 %Іо, отн. ед. 



Рис. 1. Зависимость отношения интенсивностей массовых 
линии ионов изотопов молибдена 92 Мо и Ю0 Мо от их 
энергии: 1) стандартное изотопное отношение, 2) из 
работы [6], 3) наши измерения [9, 10] 

Результаты наших измерений подтвердили сде- 
ланный ранее вывод о существовании пика или 
плато в зависимости / ІЙ (Б), однако, вывод о корре- 
ляции энергетического положения пика с массой 
элемента не подтвердился. Вместе с тем, действи- 
тельно наблюдается уменьшение максимального 
значения отношения /и^Е) к стандартному значе- 
нию отношения распространенностей при умень- 
шении массы элемента. Кроме того, для изотопов 
одного элемента, справедливо следующее утверж- 
дение: чем больше разница количества нейтронов в 
ядрах, тем больше величина ИЭ в ВИЭ. 

Из установленных к настоящему времени зако- 
номерностей наиболее надежными являются сле- 
дующие: большая вероятность ионизации атома 
легкого изотопа по сравнению с тяжелым и обрат- 
ная зависимость величины эффекта от массы в ря- 
ду элементов. Последняя закономерность, очевид- 
но, коррелирует с относительным изменением мас- 
сы и объема ядра при одинаковой величине дефи- 
цита или излишка нейтронов в ядрах при переходе 
от легких элементов к тяжелым. Данное обстоя- 
тельство указывает, что ИЭ ВИЭ связан с различи- 
ями строения электронных оболочек атомов-изо- 
топов. Различия в строении электронных оболочек 
обусловлено изотопическим смещением (ИС) 
энергетических уровней электронов в атомах. ИС в 
основном зависит от двух факторов: конечности 
массы и пространственной протяженности ядра. 

ИС электронных уровней, связанное простран- 
ственной протяженностью ядра обусловлено тем, 
что волновая функция внешнего ^-электрона пе- 


рекрывается с волновыми функциями ядра. В этом 
случае сравнительно велика вероятность нахожде- 
ния электрона в области пространства, занимаемой 
ядром (т.е. в объеме ядра) [24, 25]. Пространствен- 
ная протяженность ядер характеризуется среднек- 
вадратичным (протонным) радиусом ядра <г 2 >. Ве- 
личина протонного радиуса зависит от количества 
нейтронов в ядре (см. рис. 2). Следовательно, от 
изотопа к изотопу для одного и того же элемента 
варьируется перекрытие волновых функций внеш- 
них электронов и ядер. Таким образом, размеры 
ядер-изотопов через волновую функцию внешнего 
электрона влияют на величину 1-го потенциала ио- 
низации атома, а следовательно и на вероятность 
ионизации атома при распылении. Это влияние мо- 
жет проявляется через различия: 1) величин потен- 
циалов ионизации, 2) длины взаимодействия (у) 
атомов с поверхностью при отлете, изменяющихся 
вследствие неодинаковых величин <г 2 > для ядер- 
изотопов, поскольку у~а~<г 2 >, где а -радиус атома. 



Рис. 2. Зависимость среднеквадратичного радиуса ядра от 
числа нейтронов для изотопов рубидия [26] 


Наблюдаемые изотопические сдвиги частот 
спектральных линий некоторых оптических перехо- 
дов достигают КН... КН эВ [24]. Энергии ионизации 
(I) атомов-изотопов должны различаться на величину 
ИС, наблюдаемых в оптике. Вероятность ионизации в 
ВИЭ экспоненциально зависит от I. Это позволяет 
оценить относительную разницу вероятностей иони- 
зации атомов-изотопов. В соответствии с [26] вероят- 
ность ионизации ІГ определяется соотношением: 


К + (Е)~Е"-ех р 



( 1 ) 


где Е - кинетическая энергия иона, I - потенциал 
ионизации атома, ср - работа выхода электрона, 
с = п/2уЬ, Ѣ - постоянная Планка; п = у а 8/АукТ„ 
у а = а р /ѵ, а р - величина порядка постоянной решет- 
ки, ѵ - скорость отлетающего атома, 5 - константа, 
к - постоянная Больцмана; 7] - температура элект- 
ронной подсистемы в области развития каскада; 
\/у а - длина взаимодействия между отлетающим 
атомом и поверхностью. На основании ( 1 ), величи- 
на относительной разности вероятностей иониза- 
ции двух атомов, различающихся потенциалами 
ионизации на величину А/ и длиной взаимодей- 
ствия на величину ДА = Д(1/у„), равна: 
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К(Е) 


К -к: 

к: 


Я | _ ]?і & А*і 

к 


■ехр 


А/ -к. 


( 2 ) 


где К\ и к 2 - константы размерности (ниже при рас- 
четах полагалось к 1 =/с 2 =1). 

К. отн. ед. 



Рис. 3. Сравнение расчетных и экспериментальных зависи- 
мостей относительных вероятностей ионизации изо- 
топов вторичных ионов от их энергии: 1) эксперимен- 
тальные данные, обработанные в соответствии с (4); 
2) и 3) расчет, соответственно, по формулам (2) и (3) 


На рис. 3 представлены теоретические кривые 
Я(Е), рассчитанные по (2) - кривая 2, а также по 
формуле (кривая 3): 

Я(Е) = ^ = Е & *-\. ( 3 ) 


При расчетах полагалось, что /с , ■ А;/= /с э • А /~ К) 4 . 
Кроме того, на рис. 3 (кривая 1) представлена из- 
меренная экспериментально для изотопов серебра 
107 А§ и Ш9 А§ величина 


Я\Е) = 


А М 
М 


Я 


Я 



(4) 


где Я экс „=ГЦЕ)/Іи(Е), Іг (і=Ь,Н) - интенсивности 
массовых линий изотопов в соответствующих 
энергетических спектрах, Я ш - стандартное отно- 
шение изотопных распространенностей, АМ/М - 
относительная разность масс изотопов. Формула 
(3) получена аналогично (2) с использованием за- 
висимости Я*~Е‘ вместо (1). Такие зависимости 
также получены в ряде моделей ионизации вторич- 
ных атомов. Таким образом, кривые 2 и 3 рис. 3 ха- 
рактеризуют проявление в ИЭ ВИЭ двух процес- 
сов, один из которых описывается степенной, а 
другой экспоненциальной функцией, а кривая 1, 
возможно, характеризует суперпозицию этих про- 
цессов. На рис. 4 представлены зависимости вели- 
чин К(Е) (4) для ряда элементов, полученные в на- 
ших экспериментах и построенных по результатам 
работы [6]. Сравнение рис. 3 и рис. 4 показывает, 
что форма каждой из экспериментальных кривых в 
основных своих чертах соответствует одной из тео- 
ретических кривых. Различия в положении острых 
минимумов на экспериментальных и теоретичес- 
ких кривых можно объяснить, во-первых, экспери- 
ментальными погрешностями, основная из кото- 


рых связана с использованием энергоанализаторов 
с относительно широким энергетическим окном 
(~1 эВ), во-вторых, недостаточной обоснован- 
ностью использованных при расчетах величин 
АЯ=А( 1/т - а )=10^. Абсолютные значения сравнивае- 
мых экспериментальных и теоретических величин 
совпадают лишь в области энергий ЮЛ.. 10 2 эВ. 


К', отн. ед. 




Рис. 4. Экспериментальные зависимости (7? эш //? ст ,)— 7 от энер- 
гии вторичных ионов: а) по результатам наших изме- 
рений для изотопов: 1) 24 Мд и 26 Мд , 2) 92 Мо и 94 Мо, 
3) 50 Сг и 52 Ст, 4) " 4 5п и ІК 5 п; б) по результатам работы [6] 
для изотопов: 1) 5Л и 5І 30 , 2) О 50 и О 52 , 3) Мд 24 и Мд 26 

Проведенное сопоставление теории и экспери- 
мента, несмотря на значительные расхождения в по- 
ложении особенностей кривых по шкале энергий, 
указывает на возможность резонанса в зависимости 
№(Е) для атомов в области энергий 0...10 2 эВ. 

2. Изотопные эффекты при отражении ионных пучков 

Изотопные эффекты в потоке обратно рассеян- 
ных от поверхностей металлов ионов изучены в ос- 
новном для пар ионов 3 Не + - 4 Не + и 20 Ке + - 22 Ые + 
[11, 12]. Исследования проводились методом двой- 
ной спектроскопии ионного рассеяния (ДСИР). В 
ДСИР компоненты указанных пар присутствуют в 
первичном пучке в равных концентрациях. Исполь- 
зовался также метод поочередного напуска в разряд- 
ную камеру ионного источника указанных изото- 
пов. Энергия первичных ионов варьировалась от 
-200 до -3000 эВ, плотности тока - 0,1. ..500 нА, уг- 
лы падения и углы отражения - от нормальных до 
скользящих. На рис. 5 представлены сводные дан- 
ные из работы [12] по изотопному эффекту для пары 
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ионов Не - Не". Выход рассеяния (У) на данном 
рисунке дан в произвольных единицах, хотя относи- 
тельные разности интенсивностей пар ’Не + -Не + для 
каждого из металлов соответствуют эксперименту. 
Начала осей абсцисс для каждого из металлов 
(7п, РЬ, Іп, §п ) показаны чертой под соответствую- 
щим символом. По оси абсцисс отложены обратные 
скорости ( 1/ѵ=т/2Ѵе , т не- масса и заряд иона, 
V - ускоряющий потенциал). Хорошо развитая 
структура максимумов на кривых ДСИР, для случа- 
ев РЬ, Іп и §п на рис. 5, объяснена в [27] как резуль- 
тат конкуренции между ростом сечения рассеяния и 
вероятности нейтрализации, по мере уменьшения 
энергии пучка. В области низких значений обрат- 
ной скорости для всех комбинаций ион-мишень ток 
рассеянных ионов ’Не + выше, чем ток 4 Не + . Это 
ожидаемо, легкие бомбардирующие частицы имеют 
большее сечение рассеяния по сравнению с тяжелы- 
ми при прочих равных условиях (например, силовое 
поле, вероятность нейтрализации). При высоких 
значениях обратной скорости имеется тенденция к 
тому, что ток рассеянных ионов 4 Не + превышает ток 
рассеянных ионов 3 Не + для всех изученных систем. 
Это возможно, если при низких энергиях ионы 4 Не + 
нейтрализуются с гораздо меньшей вероятностью 
по сравнению с 3 Не + . Очевидное различие между 
этими двумя изотопами, имеющими одну и ту же 
скорость, состоит в различии их импульсов, кото- 
рое, в свою очередь, создает небольшие различия 
(намного меньше 1 ед. Бора) в их траекториях столк- 
новения. Поскольку считается, что основные про- 
цессы нейтрализации (оже-нейтрализация и резо- 
нансная туннельная нейтрализация) происходят 
сравнительно далеко от поверхности (около 4 ед. Бо- 
ра), трудно понять, как такие малые различия траек- 
торий могут влиять на относительную вероятность 
нейтрализации ионов изотопов. 



Рис. 5. Зависимость выхода упруго рассеяния ионов 3 Не + и 
4 Не 3 У от их обратной скорости для поверхностей цин- 
ка, свинца, индия и олова; угол рассеяния — 90° [12] 


Отметим также, что более тяжелые ионы, вызы- 
вают более сильное модифицирование поверхнос- 
ти в результате имплантации, поверхностного 
травления и очистки. Однако, эти процессы не да- 
ют приемлемого объяснения характера кривых 
ДСИР при низких энергиях. 

Осцилляции кривых ДСИР рис. 5 для ионов ге- 
лия, рассеиваемых на поверхности свинца, индия и 
олова, были приписаны [28] квазирезонансному 
электронному переносу. При рассмотрении этих 
кривых можно заметить, что положение максиму- 
мов и минимумов при рассеянии 3 Не + и 4 Не + на од- 
ной и той же поверхности соответствуют приблизи- 
тельно одним и тем же значениям обратной скорос- 
ти [28]. Более детальное рассмотрение показывает, 
что в действительности в области низких энергий 
наблюдаются значительные смещения в положени- 
ях максимумов и минимумов. В случае бомбарди- 
ровки поверхности свинца (рис. 5, в) изотопный 
сдвиг пика при значении обратной скорости около 
12-10 8 с/см согласуется с данными по исследованию 
угловых распределений [29]. То есть, увеличение уг- 
ла рассеяния при постоянстве энергии и массы 
бомбардирующей частицы оказывает такое же вли- 
яние, что и увеличение ее массы при постоянстве 
скорости и угла рассеяния. В обоих случаях умень- 
шается расстояние наибольшего сближения. Инте- 
ресно отметить, что сдвиги низкоэнергетического 
пика для индия и олова (рис. 5 , с, сі) имеют другой 
знак по сравнению со сдвигами для свинца. 

Из рис. 5 также видно, что высоты пиков, для 
3 Не + и 4 Не + не идентичны и не пропорциональны 
друг другу. Некоторые различия высот пиков в низ- 
коэнергетических спектрах от олова, вероятно, 
обусловлены различиями общей формы кривых. 
Наиболее отчетливо это демонстрируют неосцилля- 
торные данные для цинка (рис. 5, а). В целом кри- 
вые имеют форму параболы с синусоидальной моду- 
ляцией. Исключением из этого правила является 
кривая рассеяния ионов 3 Не + от РЬ, имеющая не- 
большой пик при значении обратной скорости око- 
ло 8,5- ІО 8 с/см, связанный, по-видимому, с квазире- 
зонансным электронным переносом. Количествен- 
но, токи отраженных ионов легких изотопов могут в 
случае 3 Не + и 4 Не + примерно в 2 раза превышать то- 
ки тяжелых изотопов. В случае 20 № + и 22 № + это пре- 
вышение составляет около 20 %. То есть, чем тяже- 
лее элемент, тем меньше величина эффекта. Это 
коррелирует с величиной относительной разности 
масс изотопов (как и в случае ИЭ в ВИЭ, см. выше). 

3. Изотопные эффекты при ионной имплантации 

Результаты наших исследований изменения изо- 
топного состава меди, имплантированной в никель, 
опубликованы в [21-23]. Ионную имплантацию ме- 
ди в поликристаллический никель со средним разме- 
ром зерна около 20 мкм проводили на имплантере 
"Диана-2" [30], оснащенном вакуумно-дуговым ион- 
ным источником, работающим в импульсно-перио- 
дическом режиме с длительностью импульса 
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150.. .200 мкс и частотой следования 50 Гц. Имплан- 
тация осуществлялась при давлении ~6- 1 0 3 и 
410 2 Па с плотностью ионного тока в импульсе 
125 мкА-см~ 2 и средней плотностью тока на мишени 
порядка 1 мкА-смЛ Дозы имплантации - 2- ІО 16 и 
2- ІО 17 ион-см~ 2 , разница во времени между набором 
указанных доз около 6 ч. Ионы ускорялись в поле 
50 кВ. Образцы никеля, предварительно электрохи- 
мически полированные, укреплялись металлическим 
держателем так, чтобы обеспечивался плотный кон- 
такт с поверхностью металлической подложки. Тем- 
пература образца в процессе имплантации не превы- 
шала 50 °С. На рис. 6 приведены профили распреде- 
ления изотопов меди 65 Си и 63 Си по глубине образца 
никеля и суммарный 65 Си+ 63 Си профиль после имп- 
лантации дозой 2- ІО 17 ионсм -2 при давлении 6- 10 3 Па. 
Видно, что характер распределения изотопов разли- 
чен: спадающая кривая для изотопа 65 Си, и нарастаю- 
щая - для изотопа 63 Си, т.е. в поверхностном слое, 
примерно до 100 нм происходит накопление изотопа 
65 Си. И только после выхода на максимум интенсив- 
ности линии изотопа 63 Си начинается выравнивание 
изотопного состава, приближаясь к их природному 
соотношению на глубине 250 нм. Для большей наг- 
лядности по данным рис. 6 на рис. 7 представлены 
отклонения концентраций (С,) изотопов меди по глу- 
бине, рассчитанные как С=//І/. (і - изотоп, I - ин- 
тенсивность), от их природных распространенностей 
(Д). Видно, что на глубине примерно 330 нм получен- 
ные значения распространенностей изотопов имп- 
лантированной меди совпадают с природными зна- 
чениями. Максимальное отклонение от природной 
распространенности изотопов (около 54 %) наблюда- 
ется на поверхности. По данным рис. 6 и 7, выявля- 
ются 3 характерные области изменения изотопного 
состава: 1 -от0 до 80.. .100 нм, 2 - от 100 до 250 нм, 3 
- глубже 250 нм. Видно, что изменение изотопного 
состава наиболее значительно в приповерхностном 
слое. Здесь наблюдается нарастание концентрации 
имплантированных ионов с выходом на максимум на 
глубине -100 нм (рис. 6). Области резкого снижения 
концентрации имплантированной меди в слое 

120.. .250 нм соответствует переход к природному со- 
отношению распространенностей изотопов. И, нако- 
нец, диффузионному "хвосту" на профиле импланти- 
рованной меди соответствует область выравнивания 
изотопного состава с природным. Таким образом, 
максимальные нарушения изотопного состава про- 
исходят на участке наибольших потерь энергии имп- 
лантированных ионов и максимального накопления 
концентрации имплантированной меди. 

На рис. 8 представлены кривые, аналогичные 
рис. 7, но для дозы 2Т0 16 ион-смТ Из сравнения 
рис. 7 и 8 видно, что характер послойного изменения 
изотопного состава с увеличением дозы не изменяет- 
ся, но для более высокой дозы изменения изотопно- 
го состава распространяются на большую глубину. 
Это, возможно, связано с увеличением общей кон- 
центрации имплантированной меди и возрастанием 
ее концентрации на большей глубине в результате 
радиационно-стимулированной диффузии. Экспе- 


римент был повторен при дозе ионной имплантации 
2- ІО 17 ионсм 2 при давлении —410 2 Па. В этом случае 
эффект трансформации изотопного состава меди ка- 
чественно не изменяется, однако, обнаруживаются 
различия в форме профилей изотопов имплантиро- 
ванной меди (рис. 8). Они становятся более узкими, 
стягиваясь к поверхности образца. 


У,отн.ед. 



Рис. 6. Профили распределения изотопов 65 Си и 63 Си по глу- 
бине образца никеля и суммарный 65 Си+ 63 Си профиль 



Рис. 7. Отклонения концентраций (0) изотопов меди по глу- 
бине от их природных распространенностей (Я,) 


С -К, % 



Рис. 8. Отклонения концентраций изотопов меди по глубине 
образца никеля от их природных распространенностей 

Объяснение трансформации изотопного соста- 
ва имплантированной меди невозможно в рамках 
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традиционных механизмов диффузии. Известно, 
что отношение коэффициентов диффузии изото- 
пов обратно пропорционально корню квадратному 
из отношения их масс, т.е. -1,02. Наблюдаемый 
эффект в несколько раз выше. 

Заключение 

Представленные в первой части статьи экспе- 
риментальные результаты касаются лишь немногих 
техногенных воздействий на поверхность, приво- 
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